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Abstract 

The geometric measure of entanglement for a symmetric pure state with nonnegative 
amplitudes has attracted much attention. On the other hand, the spectral theory of 
> ■ nonnegative tensors (hypermatrices) has been developed rapidly In this paper, we 

f^ ■ show how the spectral theory of nonnegative tensors can be applied to the study of the 



geometric measure of entanglement for a pure state with nonnegative amplitudes. Es- 
pecially, an elimination method for computing the geometric measure of entanglement 
© ; for symmetric pure multipartite qubit or qutrit states with nonnegative amplitudes is 

given. For symmetric pure multipartite qudit states with nonnegative amplitudes, a 
numerical algorithm with randomization is presented and proven to be convergent. We 
show that for the geometric measure of entanglement for pure states with nonnegative 
amplitudes, the nonsymmetric ones can be converted to the symmetric ones. 
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1 Introduction 

The quantum entanglement problem is regarded as a central problem in quantum infor- 
mation [20J, and the geometric measure is one of the most important measures of quantum 

entanglement [H0[TTl[Tll[2Tl[3Ol|3T] ^ W&S ^ rS ^ P ro P osec ^ ^y Shimony [30] and generalized 
to multipartite systems by Wei and Goldbart [31] , and has become one of the widely used 
entanglement measures for multiparticle cases [71fTT | [T ^ [25] . Among them, the study on pure 
states with nonnegative amplitudes attracted much attention recently. Wei and Goldbart 
conjectured in [3TL Section II. A] that the nearest separable state for a symmetric state can 
be chosen to be symmetric. Hayashi, Markham, Murao, Owari and Virmani [11] proved the 
conjecture for the special case of symmetric states with nonnegative amplitudes and that 
the nearest separable state can be chosen with nonnegative amplitudes in this situation, 
and Hiibener, Kleinmannn, Wei, Gonzalez- Guillen and Guhne [2] proved the conjecture 
completely. The computation of the symmetric pure states with nonnegative amplitudes 
was carried out by Wei and Goldbart [31] for some ground states, and systematically for 
symmetric pure multipartite qubit states by Chen, Xu and Zhu [7j. For general evalua- 
tions of the geometric entanglement for symmetric pure states, please see Orus, Dusuel and 
Vidal [21], Chen, Xu and Zhu [7], and references therein. 

The central problem of the computation of the geometric measure is to find the largest 
entanglement eigenvalue [12j[25j[3T]. Mathematically, the quantum eigenvalue problem is a 
generalization of the singular value problem of a complex matrix [T2JE5]. There have been 
several generalizations of singular values / eigenvalues of matrices to tensors (hypermatrices) 
recently j5j[T8l[22]. These form the spectral theory of tensors, please see [23] and references 
therein for the state of the art. In this paper, we investigate the geometric measure of pure 
states with nonnegative amplitudes through the spectral theory of nonnegative tensors. 

More specially, we establish a connection between the concept of Z-eigenvalues of tensors 
and the quantum eigenvalue problem. We show that the geometric measure of a symmetric 
pure state with nonnegative amplitudes is equal to the Z-spectral radius of the underlying 
nonnegative tensor (Theorem I3.ip . Based on this connection, a method in views of variable 
elimination [251 EZ] for computing the geometric measure of entanglement for symmetric 
pure multipartite qubit or qutrit states with nonnegative amplitudes is given. For the qubit 
case, it is an alternative to the method in [7J Section II. A]. For the qutrit case, it is new 
and gives an analytical derivation of the geometric measure of entanglement for such states. 
Here a qutrit means a system in a three dimensional Hilbert space, and a qudit used in the 
sequel is for a higher dimeniaonl case, as in 



For symmetric pure multipartite qudit states with nonnegative amplitudes cases, a nu- 



merical algorithm with randomization (Algorithm IA.2J) is presented and proven to be con 



vergent. The method is based on the shifted higher order power method (Algorithm lA.lj) 



analyzed in [To]. We show that if the initial points are randomly chosen from the intersection 
of the positive orthant and the unit sphere, then with a positive probability the algorithm 
finds the geometric measure of such states (Theorem lA.ip . The probability of the conver- 
gence of the algorithm is determined by the distance between the Z-spectral radius and the 
second largest nonnegative Z-eigenvalue of the underlying tensor. 

For nonsymmetric pure states with nonnegative amplitudes, similar results are estab- 
lished as well. Moreover, we show that for the geometric measure of entanglement for pure 
states with nonnegative amplitudes, the nonsymmetric ones can be converted to the sym- 
metric ones. Consequently, the numerical methods established above are applicable. 

The rest of this paper is organized as follows. The definitions and some basic facts of 
the geometric measure and the Z-eigenvalues of tensors are presented as preliminaries in the 
next section. In Section 3, a connection between the geometric measure of a symmetric pure 
state with nonnegative amplitudes and the theory of the Z-spectral radius of a nonnegative 
tensor is established. The computational issues are discussed. The details of the numerical 
algorithm for multipartite qudit states are put in Appendix. In Section 4, a connection be- 
tween the geometric measure of a nonsymmetric pure state with nonnegative amplitudes and 
the spectral theory of nonnegative multilinear forms is established. The paper is concluded 
with some final remarks in Section 5. 



2 Preliminaries 

In this section, some preliminaries of the geometric measure of quantum entanglement 
and the Z-eigenvalues of tensors (hypermatrices) are presented. 

2.1 Geometric measure 

An m-partite (m > 3 in this paper) pure state |\?) of a composite quantum system can be 
regarded as a normalized element in a Hilbert tensor product space % = (S^i^An where 
the dimension of Hk is dk for k = 1, . . . ,m. A separable m-partite state |<&) G % can be 
described by |$) = (g)™ =1 \<f>W) with \<f>W) G U k and |||0 (fc) )|| = 1 for k = 1, . . . ,m. A state 
is called entangled if it is not separable. In this paper, only pure states are considered. 

For a given m-partite pure state |\I/) G "H, one considers its nearest separable state 



< ^ ) ) = ®k=i lv ) m terms of the maximal overlap: 



G(V)= max |(M)|. (1) 



The geometric measure is then defined as [31] 

£ G (|*}) = 1-G(*f. 

It is shown that the maximal overlap in (prj is equal to the largest entanglement eigenvalue 

<*l (<8W fe I0 (i) >) = W\ 

®^ fc (0 (i) l) *> = A|0«>, (2) 

|||0( fe ))|| = l,fc = l,...,m. 

A state |f ) G ?i = (S^fcLi^fc is called nonnegative if there exist orthonormal bases 

{I e i fc) »f=i for ^ such that o<i--im : = (*l (k?) ' ■ • \ e t ] )) ^ ° for a11 H = !> • • -A and 

j = 1, . . . , to. The di x • • • x d m mutliway array consisting of Oj 1 ...i m is denoted by Ay. 

When "Hi = ■ ■ ■ = Tim, Ay is symmetric if and only if |\l/) is symmetric in the sense of 
quantum information [TTJ [131 [20] . The geometric measure of symmetric states attracted 
much attention recently [?t [TTl[Ti], l2~T] . 

When |\I>) is symmetric, ([T]) reduces to [H] 

G(*)= max |<¥|$)|. (3) 

2.2 Z-eigenvalues of a tensor (hypermatrix) 

For a tensor (or hypermatrix) T of order m and dimension n with m, n > 2, we mean a mul- 
tiway array consisting of numbers £i r ..j m G R for all z'j G {1, . . . , n} and j G {1, . . . , m}. The 
set of all m-th order n dimensional tensors is denoted by ~R m ' n . Given a vector x G C n , define 
Tx m_1 as an n- dimensional vector with its z-th element being V* n , £,•,„...? x,„ ■ • -x,- . 

Z-eigenvalues of tensors were introduced by Qi [22] • Suppose that T is a real tensor, i.e., 
T G R m,n . A number A G 1 is called a Z-eigenvalue of T, if it, together with a nonzero 
vector xGR", satisfies 



Tx m_1 = Ax 

(4) 
xx = 1. 

x is then called an associated Z-eigenvector of the Z-eigenvalue A, and (A, x) is called a 
Z-eigenpair. Obviously, A = Tx m := ST 1 ,...,i m =i^i-im x «i ' ' " x * m for a Z-eigenpair (A,x) 



of T. A tensor T G R m,n is called nonnegative, if ti x ...i m > for all ij G {!,..., n} and 
j G {!,..., m}. 



Many interesting results on Z-eigenvalues of tensors were obtained very recently 
[TMn?P 2 P 3] . especially, for nonnegative tensors. These results give insights on the behaviors 
of the Z-eigenvalues and powerful numerical algorithms for computing the Z-spectral radius 
of a nonnegative tensor, please see [6] and references therein. 

3 A symmetric pure state with nonnegative ampli- 
tudes 

In this section, we establish a connection between the geometric measure of entanglement 
for a symmetric pure state with nonnegative amplitudes and the Z-spectral theory of non- 
negative tensors. Based on this connection, the computation of the geometric measure for 
such states is investigated. 

3.1 A connection 

The following result was established in [TU Theorem 1]: 

Proposition 3.1 If\^f) G % is symmetric and nonnegative with the underlying orthonormal 
basis {|e;)}" =1 , then |$) = |0) 0m in ([3]) can be chosen with (ej|0) > for all i — 1, . . . ,n. 

Denote by R™ the nonnegative orthant of R™, R+ + the interior of R™ , and <S n_1 the unit 
sphere in R n . Then, we have the following result. 

Corollary 3.1 If \^) G % is symmetric and nonnegative, then 

Gm = max A^ m . (5) 

xeR^ruS 11 - 1 

Proof. Under the orthonormal basis which makes .4^ nonnegative, we have that 

GO) = max |ApX m |. (6) 

x ff x=l 

Here the superscript means conjugate transpose. For any x with x^x = 1 being an 
optimal solution for problem flB]), let x = |x| be the componentwise module of x. Then, 
x T x = 1 and 

A^ m < |Ai,x m | < ^|x| m = A^ m . 



Here the first inequality follows from the facts that x is optimal and x is feasible for (JSJ), 
and the second from the fact that .4.^ is nonnegative. Consequently, the result (j5]) follows. 

□ 

To establish a connection, we present some basic results on the Z-eigenvalues of nonneg- 
ative tensors. The following concept is important for nonnegative tensors. T = (ti 1 i 2 ...i m ) is 
called reducible if there exists a nonempty proper index subset I C {1, . . . , n} such that 

If T is not reducible, then T is called irreducible. Denote by Z{T) the set of all Z-eigenvalues 
of tensor T and q{T) : = max{|A| | A G Z(T)}. 

Proposition 3.2 Let T G K mn . 

(a) Every symmetric tensor T has at most m _ 2 ~ Z-eigenvalues. 

(b) IfT is nonnegative, then there exists a nonnegative Z-eigenpair (Ao,x^), i.e., Ao > 
and x(°) e M" D S n ~ 1 . IfT is furthermore irreducible, then Ao > and x^ G K+ + - 

(c) IfT is nonnegative and symmetric, then g(T) G Z(T) and 

g(T) = max Tx m = max Tx m 
xe-S™- 1 xeM^n^™- 1 

Proof, (a) follows from [21 Theorem 5.6], (b) follows from [6] Theorems 2.1 and 2.2], and 
(c) follows from [5J Theorem 3.10]. □ 

We call o(T) the Z-spectral radius of tensor T [5]. For nonnegative T G R m ' n , let A(T) 
be the set of A > such that A together with some x G M" fl 5 n_1 is a Z-eigenpair of 
T. Then, g{T) = max{A | A G A(T)} by Proposition 13.21 Based on Corollary 13.11 and 
Proposition 13.2} we now establish the connection. 

Theorem 3.1 // |\l/) G % is symmetric and nonnegative, then 



3.2 Computation 

Theorem 13.11 shows that the geometric measure of symmetric pure states with nonnegative 
amplitudes [HEH] can be computed through finding the Z-spectral radii of the underlying 



nonnegative tensors. In this subsection, based on Theorem 13.11 the computation of the 
geometric measure of such states is discussed. 

3.2.1. Multipartite qubit states. For symmetric pure multipartite qubit states 
with nonnegative amplitudes, [7J Section II. A] converts the geometric measure G(ty) into a 
polynomial rational fraction in one variable. By the derivatives, G(\I/) can be computed. On 
the other hand, with [261 Theorem 1], G(^) can be computed through finding the roots of 
an univariate polynomial of degree m. These two methods are somewhat different but both 
of them are the variable elimination methods. 

3.2.2. Multipartite qutrit states. We consider multipartite qutrit states in this 
part. The separability and measure of qutrit entanglement were discussed by Caves and 
Mulburn [3], and Hassan and Joag [TO] , and attracted much attention. 

Based on the variable elimination method [261 Theorem 3] and [271 Appendix], we can 
compute the geometric measure of entanglement for symmetric pure multipartite qutrit 
states with nonnegative amplitudes. 

In the sequel, we present the details. Let %\ = ■ ■ ■ = Tim, d\ = ■ ■ ■ = d m = 3, and 
"H := <^) 7 ^ =1 'Hk- Given a symmetric pure state |\1/) G "H, if |\1/) is nonnegative, i.e., there 
exists a basis {|ej)} 3 =1 such that the tensor Am is nonnegative, then G(^/) is equal to the 
Z-spectral radius of the tensor .4^ by Theorem 13.11 

Denote by Oi 1 i 2 ...i m the (i\, i 2l . . . , z m )-th element of the tensor ^4^, we have the system 
of the Z-eigenvalue equations (B| of the tensor Aqi is: 



(7) 



Note that A = *4^x 3 for any Z-eigenpair (A,x) of ^4*. We now give the following 
algorithm. 

Algorithm 3.1 (Symmetric multipartite qutrit states with nonnegative amplitudes) 

Step Input data A® with elements a ili2 ,,, im . Set 11(^4^) as the empty set. 

Step 1 If O21...1 = ct3i...i = 0, then A = an...i is a Z-eigenvalue with a Z-eigenvector x : = 
(1,0,0). In this case, put A = Ou,„i into 11(^4^). 



■^12, ■ 


.,j m = l a lJ2---«m X *2 ' 


X *m 


= Axi, 


E 3 


.,j m = l a 2J2...«m X *2 ' 




= Ax 2 , 


E 3 


.,i m = l a 3J 2 ...«m X *2 ' 


■ -X,' 


= Ax 3 , 


> 


x l "r x 2 "r x 3 




= 1. 





.,i m =l a l«2---*m X *2 ' 




= Axi, 


E 2 


.,i m =l a 2J2...Jm X «2 ' 


X *m 


= Ax 2 , 


E 2 


.,j m = l a 3i2...«m X *2 ' 




= o, 


y 


X 2 + X 2 




= 1. 



Step 2 In this step, we consider Z-eigenvectors with X3 = but x 2 7^ 0. If X3 = and 
x 2 7^ 0, then (JTJ) becomes: 



(8) 



(i) Multiply the first equation of OH]) by x 2 and the second by x l5 then we get 

2 2 

x 2 / j aij 2 ...j m Xj 2 • • -Xj m = Xi / j a 2 j 2 ...j m Xj 2 • • -Xj m . 

i2i-.'m = l i2,-.'m = l 

Divide it by x™, the third equation of (JSJ) by x™ -1 , and set £ = ^. Then, we get 
two polynomial equations in variable t as f(t) = g(t) and h(t) = 0. 

(ii) For all nonnegative t such that both f(t) — g(t) = and h(t) = 0, we have 
x := f TfT^ i /THf ; j is a Z-eigenvector of *4^. The corresponding Z-eigenvalue 
is A = *4*x 3 . Put these A into 11(^4^). 

Step 3 In this step, we consider Z-eigenvectors with x 3 7^ 0. If x 3 7^ 0, then we have the 
following: 

(i) Multiply the first equation of (JTJ) by x 3 and the third by x l5 then we get 

3 

i2---im X «2 ' ' ' X im X l / j ^3l2 —im X «2 ' ' ' X «m • 

i2,...,im=l 

Similarly, we have 

3 3 
X 3 / G 2i2...im X i2 ' ' ' X im = X 2 / _ , a 3i2---im X i2 ■ ■ ' X im- 

J2v:*m = l I3,...,i m =l 

Divide them by x™ respectively, and set u = — and v = — . Then, we get two 
polynomial equations in variables u and v as f(u,v) = and g(u,v) = 0. 

(ii) Write them in univariate polynomial equations in the variable u with coefficients 
as univariate polynomials in the variable v as 

f(u,v) = a (v)u m + ai^u" 1 " 1 H ha m (v)=0 } 

g(u,v) = b {v)u m - l + b l {v)u m - 2 + --- + b m - l {v) = $. 



3 

E 

42, ■■;im= 


a i 
1 



(iii) Form the (2m — 1) x (2m — 1) polynomial matrix (the Sylvester matrix [8]) 
the variable v as: 



m 



f a (v) a%(v) ■■■ a m (v) 

a (v) ai(v) ■■■ a m (v) 



M(v) :-- 





bo(v) bi(v) 
b (v) b x (v) 




b m ~i{v) 



a (v) oi(v) 


b m -i(v) 






a m (v) 





\ 



\ 











b Q (v) bx(v) ■■■ 6 m _i(u) / 



(iv) Compute the determinant of the matrix M(v) as a univariate polynomial in the 
variable v. Denote it by d(v). 

(v) For any nonnegative v such that d(v) = and then any nonnegative u such 
that (u,v) being a solution for ®, x := ( . % — =, , ",, — =, -= — L — 9 ) is a Z- 
eigenvector of .4.^. The corresponding Z-eigenvalue is A = *4^x 3 . Put all such A 
into ntyW). 



Theorem 3.2 // |\P) E 7i = (£)™ =1 'Hk is a symmetric nonnegative qutrit state and II (.4.*) 
is generated by Algorithm \3.1[ then 11(^4^) = A(Ay), and 

G($) = g{Ay) = max{A | A € II(yW)}. 



Proof. By Steps 1-3 of Algorithm 13.11 and the Sylvester theorem [8] used in Step 3, all 
the nonnegative Z-eigenvectors are considered. Consequently, II (.4^) = A(Ay). Now, by 
Proposition 13.21 and Theorem 13.11 the result follows. □ 

This result is new for computing G(*&) in the qutrit case. 

3.2.3. Z-spectra of multipartite qubit and qutrit states. In the following, we 
show the Z-spectra for GHZ, W, and inverted- W states, and a qutrit state. As in [31], 
define 



\S(m,k)):-- 



k\(m — k)\ 



nil 



Y, |r(£— (J4— i)>» 



ree r 



m—k 



here & m is the symmetric group on m elements. Examples 13.11 13.21 and 13.31 are computed 
through Theorem 1 in [26J. 



Example 3.1 The mGHZ state is defined as: 

\mGHZ) := (\S(m, 0)) + \S(m, m))) /V2. 
Under the basis {|0), |1)}, we have A m Guz £ M m ' 2 and the Z-eigenpairs are: 



i=,(l,0)), (^,(0,1)) 



and 



i ill 



and five more when m is odd: 



GM- 1 . )).^".- 1 )).^. (-*.*)) 



^T"; 775; — 775 5 an d 



1 1 



We have, G{mGHZ) = g{A m GHz) = 75 which agrees with that in [311 Section II. A], and 
it can be attained with nonnegative Z-eigenvectors which agrees with Proposition 13.11 The 
corresponding nearest separable state is |$) = |0)® m with \<p) := |0) or |1). 

Example 3.2 In this example, W state for a 3-partite qubit setting is considered. The W 
state is defined as: 

\W) := \S(3, 2)> = (|001) + |010) + |100» /V3. 

Under the basis {|0), |1)}, we have Aw £ K 3 ' 2 and the Z-eigenpairs are: 



(0,(0,1)) 



2 / /2 1 \ \ 2 /_ 2 1 

' V 3' I V 3' V 3 / / ' I 3' I V 3' V 3 



iU/i-VI M nd -I» -VW§ 



Again, G(Aw) = q(Aw) = § which agrees with that in [3TJ Section II. A], and it can be 
attained with nonnegative Z-eigenvectors which agrees with Proposition 13.11 The corre- 
sponding nearest separable state is |$) = |0)® 3 with \4>) := J ||0) + -i/||l). 



Example 3.3 In this example, inverted-iy state for a 3-partite qubit setting is considered. 
It is defined as 

\W) := |S(3, 1)) = (|110) + |101) + |011>) /VS. 

Similarly, we have A^ G M 3 ' 2 . After switching xi and x 2 , the Z-eigenvalues equations (jlj) 
of W becomes that for W. Consequently, G(W) = g(A^) = | by Example 13.21 with the 
corresponding nearest separable state being |<I>) = |0)® 3 with |0) := J ||0) + w||l). 



10 



Example 3.4 In this example, a general GHZ state [101 Eq. (9)] for a 3-partite qutrit 
setting is considered. It is defined as 

|tf) := a|lll) + /3|222) + 7 |333), a 2 + /3 2 + 7 2 = 1. 

Here {|1), |2), |3)} is the basis for each qutrit. We see that Ay G M 3,3 is nonnegative and 
symmetric when a, (3, 7 > 0. In this situation, the Z-eigenvalue equations (J7j) become: 

ax] = Axi, /3x 2 = Ax 2 , 7X3 = Ax 3 , and x 2 + x 2 + x 2 = 1. 

When a/37 = 0, the Z-spectra can be computed through Theorem 1 in [2B]. We mainly 
consider the nondegenerate case when a/37 > 0. By Algorithm l3.lt we can compute all the 
nonnegative Z-eigenpairs as: 



(a, (1,0,0)), (/3, (0,1,0)), (7,(0,0,1)), (^g£=, (-^ ^=, 0) j , 

" 7 , ( , 7 , 0, ° ) 1 , ( jtL , (0, , 7 , / ) ) , and 

(Sf 1 , (?, f , f )) with r := v/a2/32 + /3V + « 2 7 2 - 

We see that G(\l/) = £>(-A*) = max{a,/3,7}. The corresponding nearest separable state is 
|$) = |</>)®3 with |0^ := |^ when q^ = a> | 2 ) w hen G{9) = (3 and |3) when G(tf) = 7. 

3.2.4. Multipartite qudit states. For symmetric pure multipartite qudit states, there 
is no analogue result of Proposition 13.21 (a) for the quantum eigenvalue problem (J2J). Conse- 
quently, the computation for the largest quantum eigenvalues of a state is very complicated 
in general. 

Nonetheless, for symmetric pure states |\l/) with nonnegative irreducible Ay, Proposition 
13.21 (a) helps to prove that the shifted higher order power method [15] is locally convergent, 
i.e., with the initial point sufficiently close to a Z-eigenvector of the Z-spectral radius, the 
shifted higher order power method converges to the Z-spectral radius (Lemma lA.ll) . Conse- 
quently, if we randomly choose initial points in R+ + n<S n_1 , then with a positive probability 
we can find the Z-spectral radius for an irreducible nonnegative tensor (Theorem lA.il) . This, 
together with Theorem 13.11 implies that Algorithm IA.2I can find the geometric measure of 
symmetric pure states |\l/) with nonnegative irreducible Ay. We present the details in Ap- 
pendix. When Ay is reducible, the technique in [13] is applicable. 



Besides the shifted higher order power method, a generalized Newton method like that 
in [17] may be developed for finding the Z-spectral radius. 



11 



4 A nonsymmetric pure state with nonnegative ampli- 
tudes 

In this section, we extend the results in the last section to nonsymmetric pure states 
with nonnegative amplitudes. To this end, we need the spectral theory for multilinear 
forms [HJUHIEB]- We first establish analogue results for nonnegative multilinear forms. 

Let A = (cLi 1 ...i m ) be a d\ x • • • x d m real tensor (hypermatrix) . a G IR is called a singular 
value of A, if it, together with x^ G M. dl n S dl ~\ . . ., x^ m ) G R dm n S dm -\ satisfies 

J^ a% 1 ...i m yil ' ' ' x £ = ax 4 > Vi fc = 1, . . . , 4, Vfc = 1, . . . , m. (10) 

The vector x( fc ) is called the mode-/c singular vector corresponding to the singular value 
a [H1I28]. Denote the largest singular value of A by a (A). 

Proposition 4.1 Let A = {a>i 1 ...i m ) be a d\ x • • • x d m real tensor. Then, 

*M) = m d ^ax , d i A^ ■ ■ ■ xW := E ■ ■ • E "---^ ' ' ' *£? '■ (") 

x( 1 )e<s d i- 1 ,...,x( m )e<s d ™- 1 z — ' z — ' 

Moreover, if A is nonnegative, then the mode-k singular vector corresponding to o~(A) can 
be chosen nonnegative. 

Proof. We see firstly from ([TO]) that for any singular value a of A with singular vectors 
x^ 1 ), . . . ,x (m ), we have a = Ax.^ ■ ■ -x^ m \ Secondly, by optimization theory, the singular 
vectors are exactly the critical points of the maximization problem (ITT]) . Hence, (fTTj) follows. 

Now, we show the second result. Suppose now that A is nonnegative. We have 

max ^x (1) • • • x (m) 

x( 1 )eK+ 1 n5 d i- 1 ,...,x( m )eR+ m n5 d "i- 1 

< max ^x (1) • • ■ x (m) 

xC 1 )s5 dl_1 ,...,x( m )e5 dm_1 

< max ^lx (1) ---x (m) . 
x( 1 )eR* L n5 d i- 1 ,...,x( m )6iR+ m n5 d "i- 1 

Here the second inequality follows from the fact that A is nonnegative. 

Consequently, 

a(A) = max ^x (1) • • • x (m) = max ^x (1) • • • x (m) . 

x ( 1 )e5 d i- 1 ,...,xf m )s5 d m- 1 x( 1 )eK^ 1 n5 d i- 1 ,...,x( m )eR dm n5 d m- 1 

12 



Hence, the optimal value o~(A) of (ITT]) can be achieved with nonnegative x^'s. Then, there 
exist nonnegative x^'s that are critical points of (ITT]) . By the correspondence of the critical 
points of ( ITT]) and the singular vectors of A, the result follows. □ 

We now generalize Corollary 13.11 to a nonsymmetric pure state with nonnegative ampli- 
tudes. 

Proposition 4.2 If\^) E H is nonnegative with the underlying orthonormal bases {\e\ }}j=i 
for k = 1, . . . ,m, then |$) = <S>™ =1 \4>^} in (JTJ) can be chosen with (e\ \<p^) > for all 
i = 1, . . . , dk and k — 1, . . . , m. Consequently, 

G(^)= max ^x (1) ---x (m) . 

x( 1 )e]R^ 1 n5 d i- 1 ,...,x( m )eiRj n n5 d «i- 1 

Proof. It is similar to that for Corollary 13. II □ 

By Propositions 14.11 and 14.21 we have the following theorem. 

Theorem 4.1 If \^) E H is nonnegative with the underlying orthonormal bases {\e\ )}jX\ 
for k = 1, . . . ,m, then 

G(y)=a(Av). 

Theorem 14.11 gives a connection between the geometric measure of nonsymmetric pure 
states with nonnegative amplitudes and the spectral theory of nonnegative tensors. 

Proposition 14.11 generalizes Proposition 13.21 (b) and (c) to the context of nonnegative 
multilinear forms. As it can be seen from the last section that Proposition 13.21 (a) plays an 
important role in the computational issues, we establish an analogue result for multilinear 
forms in the following. To this end, symmetric embedding introduced in [28] is needed. 

Let A = {di l ...i m ) be a d\ x • • • x d m real tensor and 1S4 be the symmetric embedding 
of tensor A [28, Section 2.2 ]. 1S4 is an m-th order XlfcLi ^fc dimensional symmetric tensor. 
Then, we have the following result. 

Proposition 4.3 Let A = (cii 1 ...i m ) be a &\ x • • • x d m real tensor. Then, a is a nonzero 
singular value of A if and only if T" m o is a nonzero Z-eigenvalue of Sj,. 

Proof. The "only if part follows from [281 Theorem 4.7]. 

We show the "if part in the following. Now, suppose that y := (y^ , . . . ,y^ ) T G 
M. dl+ '" +dm fl S dl+ '" +dm ^ 1 with y( fc ) e M. dk for each A; is a Z-eigenvector of S4 corresponding 

13 



to Z-eigenvalue A^O. Suppose, without loss of generality, that y^ 1 ' ^ 0. By the definition 
of 1S4, we have 



A(y (1) fy (1) 


= £y<:> 


diH hdm 

/ j \'~ > A)hi2---im~yi2 ' ' ' Yim, 






= (m-l)\Ay {1) ---y {m) 






dk 


7 j \SX)i 1 i 2 ...i rn y 

l<ij<di-\ \-dm, j^k 


'1 «/ ^m 




= A(y( fc )fy« 





for all /c = 2, . . . , m. Consequently, (y^) T y^ k ' = — for all k = 1, . . . , m. Moreover, 



diA \-d m 

/ j {SA)hi 2 ---i m yi2 ■ ■ ■ y« 

12,. ..,im = l 



[m 



d2 d m 

£2=1 im=l 



( 2 ). .„("») 



,W 



Ay^, Vzi = l,...,di. 



Let x( fc ) := y/my^ for all A; = 1, . . . , m. We then have 



(ra-1)!- 



-. C?2 dm 



X ( 2 )... X M = A 

«2...«m^l 2 ^l m A 



12=1 im=l 



m 



x !!\ v ^i = l,...,rfi. 



(1) (m) _ ■. 

fl il...im X ii ' ' ' X i m ~~ A 



1 



X; 



(k) 



\fm 
Similarly, we have 

l<ij<dj, j^tk 

This, together with f JTQj) . implies that , m A is a nonzero singular value of A. The proof is 
complete. □ 



Vu = 1 dk, V7c = 2 m. 



//? 



L, . . . , l* ft , 



Corollary 4.1 Lei „4 = (a^...^) be adiX- ■ -xd m real tensor. Then, it has at most 
singular values, here N := J^fcLi dk- 



{m-l) N -l 
m-2 



Proof. From the definition of 1S4, it is easy to see that if is a singular value of tensor 
A with singular vectors x( 1 ',...,x' m ', then is a Z-eigenvalue of 1S4 with Z-eigenvector 

( X (D) T (x(™)) T " ' 

y/m ' ■ ■ ■ ' yfi^l 

Now, the result follows from Proposition 13.21 (a) and Proposition 14.31 immediately. □ 

Corollary 4.2 If\^f) EH is nonnegative with the underlying orthonormal bases {\e\ )}^=i 
for k = 1, ... ,m, then 



G(V) = a{Av) 



\frrtf' 



ml 



■q(Sa 9 )- 
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Proof. It follows from Propositions 13.21 and I4.3[ and Theorem 14.11 immediately. □ 

So, for the geometric measure of entanglement for pure states with nonnegative ampli- 
tudes, the nonsymmetric ones can be converted to the symmetric ones. Consequently, the 
numerical methods in the last section are applicable. 



5 Conclusion 

We have established a connection between the geometric measure of entanglement for 
a pure state with nonnegative amplitudes and the spectral theory of nonnegative tensors. 
Especially, we have shown that the geometric measure of entanglement of a symmetric pure 
state with nonnegative amplitudes is equal to the Z-spectral radius of the underlying nonneg- 
ative symmetric tensor, and the geometric measure of entanglement of a nonsymmetric pure 
state with nonnegative amplitudes is equal to the largest singular value of the underlying 
nonnegative tensor. Based on this connection, the computation of the geometric measure 
is investigated in details. An analytical derivation for the geometric measure of symmetric 
pure multipartite qutrit states is given. For multipartite qudit states, an algorithm with 
randomization is proposed and its convergence is proved. We have proven that it is conver- 
gent. Independently, results in Section 4 for nonnegative multilinear forms have their own 
significance. 

Many established results and algorithms for eigenvalues of nonnegative tensors, e.g. the 
generalized Perron- Frobenius theorem [UEUHlElHlES], can De applied to the study and com- 
putation of the geometric measure of a pure state with nonnegative amplitudes. Also, there 
is a possibility that new ideas will emerge from the intersection of the geometric measure 
theory of quantum entanglement and the spectral theory of nonnegative tensors. 
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A Appendix: An algorithm for finding the Z-spectral 
radius of a nonnegative irreducible symmetric ten- 
sor 

In this section, we propose an algorithm for finding the Z-spectral radius of a nonnegative 



irreducible symmetric tensor, and prove its convergence. By Theorem 13.11 a method for 
computing the geometric measure for symmetric pure states with nonnegative amplitudes is 
then given. 

A.l The power method 

The following is the shifted higher order power method investigated in p~5] . 

Algorithm A.l (Shifted Higher Order Power Method (SHOPM)) 

Step Initialization: choose x^ G R™ + n5 n " ! and a > 0. Let k := and A := ^(x(°)) m . 
Step 1 Compute 

x (*+i) . = ^( x (fc))m-i + ax (fc) 5 x ( fc+ i) . = _g±±> an d \ k+1 ■= A(^ k+1 ^) m - 

Step 2 If A^6 k+1) ) m - 1 = A fc+1 x( fc+1 \ stop. Otherwise, let k := k + 1, go to Step 1. 

By (EH Theorem 4.4], Proposition 13.21 and [14"] . we have the following result. 

Proposition A.l Let A G W m,n be nonnegative and symmetric. For a > (m — l)g(A), the 
iterates (A&, x^) generated by Algorithm \A.l\ satisfy the following properties. 

(a) The sequence A& is nondecreasing and converges to a Z-eigenvalue A* > 0. 



(b) If there are only finitely many Z- eigenvectors, the sequence {x^} converges to a Z- 
eigenvector of A corresponding to A*. 

(c) Sequence {x^} has an accumulation point and every such point is a Z-eigenvector of 
A corresponding to A*. 

Based on Algorithm IA.ll we introduce an algorithm for finding the Z-spectral radius of 
a nonnegative irreducible symmetric tensor. 

Algorithm A. 2 (An algorithm for the Z-spectral radius of a nonnegative irreducible sym- 
metric tensor) 

Step Let k := and compute the Z-eigenvalue Ao of A through Algorithm \A.1\ Let 

Step 1 Choose N initial points in K_ + fl <S n ~ 1 . For the i-th initial point, compute the 
Z-eigenvalue fii of A through Algorithm ] A. 1\ Let fi := maxi<j<Ar /ij. 



Step 2 If fi = 1, then A := Ylo<j<k^r The algorithm is terminated. Otherwise, let \k+i 
A 



fj,, A := -t^—, and k :— k + 1, go to Step 1. 



Remark A.l By Proposition \3.2\ if A is irreducible and A G A (.4.), then A > 0. This, 
together with Proposition ] A. 11 implies that Algorithm ] A. 2\ is well defined. 

To establish the convergent of Algorithm IA. 21 we prove the following lemma first. 

Lemma A.l Let A G M. m ' n be nonnegative, irreducible and symmetric. For g{A), if x* is 
a corresponding Z-eigenvector, then there exists e > such that for any x^ ^ G K+ + fl {x G 
iS n_1 | ||x — x*|| < e} ; the sequence {A^} generated by Algorithm ] A.l] with a > (m — l)g(A) 
converges to g(A). 

Proof. If A(A) is the singleton {^(^4)}, then the result follows from Proposition IA.1I with 
arbitrary e > 0. Now, suppose that the cardinality of A (.A) is larger than one. 

Denote by X 2 (A) := max{A | A G A(A) \ {q{A)}}, and « := eW-W) _ By p ropos i t ion 



13. 2\ we see that the open set {/3 G R | \/3 — g(A) \ < k} disjoints with the union of the finitely 
many open sets {/3 G R | \fi - A| < k} for A G A(A) \ {g(A)}. Since A := ^L(x(°)) m and 
g(A) := A(-x*) m , we can choose e > such that |A - g(A)\ < k for any x (0) G R™ + fl {x G 
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S n x I ||x — x*|| < e}. Consequently, this, together with Proposition 13.2( c). implies that 
q(A) -k < A < g(A). 

By Proposition lA.lt {A^} is nondecreasing and converges to a Z-eigenvalue A* of A. As 
we can see, the only possibility is that A* = g(A). The proof is complete. □ 

Assumption A.l If the N initial points in Algorithm \A.2\ are chosen such that there is at 
least one point in the set {x G M" + niS™ -1 | ||x — x*|| < e} with e being determined by Lemma 
\A.1\ then we call such a set of initial points satisfies absolutely convergent condition 

(ACC for short). 

Theorem A.l Let A G M. m ' n be nonnegative, irreducible and symmetric. For a > (m — 
l)g(A), if the set of N initial points satisfies ACC, then Algorithm ] A. ^1 is terminated with 
k = 1 and A = g(A). In general, the sequence {Xk} generated by Algorithm \A.2\ converges 
to a positive Z-eigenvalue of A. 

Proof. The results follow from Proposition lA.lt Lemma [A. II and Assumption IA. II □ 

So, if we uniformly randomly choose initial points in M™ + fl iS n_1 , then with positive 
probability Algorithm IA.2I finds the Z-spectral radius. By Theorem IA.1I and Assumption 
IA.lt we see that the probability that Algorithm IA.2I converges to the Z-spectral radius 
under the uniformly random framework is determined by e in Lemma IA.1I In the following 
subsection, we show that e in Lemma IA.1I can be determined through k in Lemma IA.1I 
explicitly. 



A. 2 On the ACC assumption 

We in this subsection investigate the e in Lemma IA.1I It plays an important role in As- 
sumption IA.lt an d thus in Theorem IA.1I We show that there is an explicit formula for e 
in Lemma [A. II based on k := 2 ^ or an y ^ = ( a h---i m ) £ K m ' n , its Frobenius norm 

||^|| F is defined as ||.A|| F := ^Er i ,...,i m =i<...i m - 

Proposition A. 2 Let A E M m,n . For any x, y G S n ~ l , we have 

||^x m - fc -^y m - fc || F < (m-fc)||^|| F ||x-y||, V/c = 0, 1, . . . ,m. (12) 

Proof. We first show that 

H^ x (i) . . . x (m-fe) || f < ||^|| F) vjfc = 0, 1, . . . , m - 1. (13) 

Here, xW,...,^™^ G S"" 1 . 
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We show how it works for k = m — 1 in details. In this case, Axf- 1 ' G IR m 1,n . 



n / n \ 2 

ii,...,i m _l=l \j'=l / 

n n 

— / ^ / j a h im-lj 

il,...,i m _l=l j=l 

- WAWl 
Here the inequality follows from the fact that x^ 1 ) G <S n_1 . 
• In general, 

||^ x (l)... x (m-fe)|| F = p x (l)... x (™-^l) x (™-fe)|| F 

< pxW-x( m - H '| F <-< ||«4x (1) ||f 

< Mil* 

for A; = 0, . . . , m — 2. 

Now, we show the proof for (Fl2|) when m = 3. 

Mx-^ly|| F = ||^(x-y)|| 1? =||P(x-y)|| 

< l|P||||(x-y)||<||P|W|(x-y)|| 

Wl(x-y)||, 



here P is an n 2 x n matrix such that Px = vec(Ax) for all x G M n . Then, for m > 3, we 
have 

||^x m - fc -^y m - fc || F 
= ||^x m ~ fc - ^yx" 1 -^ 1 + iyx" 1 -^ 1 - ^y 2 x m - fc - 2 + • • • + Ay m - k - 2 x - Ay rn - k \\ 

< (WA^-'-'y + \\Ax m - k - 2 y\\ F + ■■■+ \\Ay m - k - l \\ F ) ||(x - y)|| 

< (m-k) \\A\\ F \\ (x - y) ||, Vk = 0, 1, . . . , m - 1. 

Here, the last inequality follows from ffT3|) . The proof is complete. □ 

Corollary A.l Let A G M. m ' n be nonnegative, irreducible and symmetric. For any x^ ^ G 
Kl+n {xG S™ -1 | ||x-x*|| < -on-}, we have |A - g(A)\ < K. 

++ l i ii ii 'm\\A\\F J ' y '' 

Remark A. 2 By Theorem \A.l\ and Corollary \ A. 11 we see that k determines the probability 
that Algorithm ] A. 2 converges to g(A) under the random framework. Note that \\A^\\f — 1 
for pure states. Then, we can analyze the probability through estimations of k for computing 
the geometric measure of entanglement. 
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